Linear magnetoresistance on the topological surface 
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A positive, non-saturating and dominantly linear magnetoresistance is demonstrated to occur in 
the surface state of a topological insulator having a wavevector-linear energy dispersion together 
with a finite positive Zeeman energy splitting. This linear magnetoresistance shows up within quite 
wide magnetic-field range in a spatially homogenous system of high carrier density and low mobility 
in which the conduction electrons are in extended states and spread over many smeared Landau 
levels, and is robust against increasing temperature, in agreement with recent experimental findings 
in Bi2Se3 nanoribbons. 

PACS numbers: 75.47.-m, 73.20.At, 73.25.+i 



I. INTRODUCTION 

It is well known that the classical magnetoresistance 
(MR) in metals or semiconductors with a closed free elec- 
tron Fermi surface increases quadratically with increasing 
magnetic field B for /iB -C 1 and saturates when fiB > 1. 
Here /t is the zero-magnetic-ficld mobility. Hence, the ex- 
traordinarily high and linear MR (LMR), which breaks 
this familiar rule, has been gaining much attention as 
soon as its discovery. In the past decade, this unexpected 
LMR has been reported in silver chalcogcnider' indium 
antimonidcf 2 - silicon^ MnAs-GaAs composite material^ 
and graphene£ 

Kapitza's linear law£ indicates that the metal shows a 
magnetoresistance linear in perpendicular magnetic field 
when it has an open Fermi surface and a mean free path 
longer than the electronic Larmor radius. Recently, an- 
other two models, irrespective of the open Fermi sur- 
face, have been constructed to provide possible mecha- 
nisms for the LMR phenomenon. Abrikosov suggested a 
quantum-limit origin of LMR for the homogenous system 
with a gapless linear energy spectrumi^ His model re- 
quires that Landau levels are well formed and the carrier 
concentration is small that all electrons occupy only the 
lowest Landau band. Alternatively, Parish and Little- 
wood developed a classical model without involving lin- 
ear spectrum^ Ignoring the concrete microscopic mech- 
anism, they attributed this unusual MR to the mobility 
fluctuations in a strongly inhomogenous system. 

Topological insulators^— (TIs) are novel materials 
with a full energy gap in bulk, while there are gapless sur- 
face states. Due to its unique band structure with only 
one helical Dirac cone and linear energy dispersion^ - — 
the surface states of the TI Bi2Se3 become an excellent 
platform for the study of quantum-limit LMR. The re- 
cent experiment in this flat surface system, however, re- 
ported that a large positive MR, which becomes very 
linear above a characteristic field of 1^2 T, was observed 
even in an opposite situation where the carrier sheet den- 
sity is high that electrons occupy more than one Landau 
levelsJ£ Moreover, they found that raising temperature 



to room temperature almost has no influence on the ob- 
served LMR. It is striking that this observation is in con- 
flict with Abrikosov's model and also with the classical 
Parish-Littlewood model. So far a reliable theoretical 
scheme capable of explaining this novel experiment has 
still been lacking. 

In this paper, we generalize the balance-equation 
approach!! to a system modeling the surface states of a 
three-dimensional TI to investigate the two-dimensional 
magnetotransport in it. We find that a positive, non- 
saturating and dominantly linear magnetoresistance can 
appear within quite wide magnetic-field range in the TI 
surface state having a positive and finite effective g- 
factor. This linear magnetoresistance shows up in the 
system of high carrier concentration and low mobility 
when electrons are in extended states and spread over 
many smeared Landau levels, and persists up to room 
temperature, providing a possible mechanism for the re- 
cently observed linear magnetoresistance in topological 
insulator Bi2Sc3 nanoribbons J£ 



II. BALANCE-EQUATION FORMULATION 
FOR MAGNETORESISTIVITY 



We consider the surface state of a Bi2Se3-type large 
bulk gap TI in the x-y plane under the influence of a 
uniform magnetic field B applied along the z direction^ 
Following the experimental observation^ we assume 
that the Fermi energy locates in the gap of the bulk 
band and above the Dirac point, i.e. the surface car- 
riers are electrons. Further, the separations of the Fermi 
energy from the bottom of bulk band and Dirac point are 
much larger than the highest temperature (300 K) con- 
sidered in this work. Hence, the contribution from the 
bulk band to the magnetotransport is negligible. These 
electrons, scattered by randomly distributed impurities 
and by phonons, are driven by a uniform in-plane elec- 
tric field E = (E x ,E y ) in the topological surface. The 
Hamiltonian of this many-electron and phonon system 
consists of an electron part rl e , a phonon part "H p h, and 
electron-impurity and electron-phonon interactions H c i 
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and "Hop: 

H = H c + H e i + H ep + H p h- (1) 

Here, the electron Hamiltonian is taken in the form 



H c — 



i 



1 



E 



in which 7T, 



(2) 

+eA( rj ) = (ttJ, ttJ), r 3 = (xj,yj), 
(Pjx,Pjy) and tTj = (crj,crj,cr|), stand, respec- 
tively, for the canonical momentum, coordinate, momen- 
tum and spin operators of the jth electron having charge 
— e, A(r) — (—By,0) is the vector potential of the per- 
pendicular magnetic held B = Bz in the Landau gauge, 
vp is the Fermi velocity, g z is the effective g-factor of the 
surface electron, and (ib = e/2mo is the Bohr magne- 
ton with mo the free electron mass. The sum index j in 
Eq. ([3J goes over all electrons of total number N in the 
surface state of unit area. 

In the frame work of balance equation approach^ - — 
the two-dimensional center-of-mass (cm.) momentum 
and coordinate P = J2jPj an d R = ^'~ 1 Y2j r 3> an d 
the relative-electron momenta and coordinates p'j = 
Pj — P/N and r'j = rj — R are introduced to write the 
Hamiltonian H c into the sum of a single-particle cm. 
part H cm and a many-particle relative-electron part H er : 

H = H cm + Her, With 



H C m = v ¥ {n x a l i - n y a x c ) + NeE ■ R, 



H 



E 
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(3) 
(4) 



In this, 77 



NeA(R) = (n x ,n y ) is the canonical 



momentum of the center-of-mass and 7r'- = p'j+eA(r'j) 

(tt^ x , TTj 1 ) is the canonical momentum for the jth relative 
electron. Here we have also introduced cm. spin op- 
erators a% = A 'V^: and a y c = AT" 1 £V <r?. The 
commutation relations between the cm. spin opera- 
tors a* and cr^ and the spin operators <rj, crj and er| 

of the jth electron are of order of 1/N: [(Tj,<jP 2 ] = 

N- 1 2ie/3 1 f3 2 p 3 a? 3 with fa, fa, fa = (a?,y,z). Therefore, 
for a macroscopic large AT system, the cm. part H cm 
actually commutes with the relative-electron part H er in 
the Hamiltonian, i.e. the cm. motion and the relative 
motion of electrons are truly separated from each other. 
The couplings between the two emerge only through 
the electron-impurity and electron-phonon interactions. 
Furthermore, the electric field E shows up only in H cm - 
And, in view of [r' ia ,p' j/3 ] = i<5«/j(% - 1/N) ~ id a p5ij, 
i.e. the relative-electron momenta and coordinates can 
be treated as canonical conjugate variables, the relative- 
motion part Her is just the Hamiltonian of A" electrons 
in the surface state of TI in the magnetic field without 
the presence of the electric field. 

In terms of the cm. coordinate R and the relative 



impurity and electron-phonon interactions can be writ- 
ten asi&12 



(5) 



Hep = M (Q> A ) <t>Q\j qR p q - (6) 

Here U(q) and M(Q, A) are respectively the impurity 
potential (an impurity at randomly distributed position 
r a ) and electron-phonon coupling matrix element in the 
plane-wave representation, and 4>q\ = bq X + &1q A with 
Bq x and bq X being the creation and annihilation opera- 
tors for a phonon of wavevector Q = (q, q z ) in branch A 
having frequency Oq\. 

The cm. velocity (operator) V is the time variation 
of its coordinate: V = R = — i[R, H] = Uf(o"c i — a% j). 
To derive a force-balance equation for steady state trans- 
port we consider the Heisenberg equation for the rate of 
change of the cm. canonical momentum 77: 



n 



-i [IT, H] = -Ne(V x B) — NeE + F + F p , (7) 



in which the frictional forces 71 and F p share the same 
expressions as given in Ref.[l9l. 

The statistical average of the operator equation (J7J) can 
be determined to linear order in the electron-impurity 
and electron-phonon interactions H i and H cp with the 
initial density matrix po = Z~ 1 e~( H * ,b+ ' H <" : " T at temper- 
ature T when the in-plane electric field E is not strong. 
For steady-transport states we have (17) = 0, leading to 
a force-balance equation of the form 







-Nev x B - NeE + /; + f p . 



(8) 



Here v = (V), the statistically averaged velocity of the 
moving center-of-mass, is identified as the average rate 
of change of its position, i.e. the drift velocity of the 
electron system driven by the electric field E, and fi and 
fp are frictional forces experienced by the center-of-mass 
due to impurity and phonon scatterings: 



f l =J2\U(q)\ 2 qn 2 (q,u J o), 
q 

f P =J2 \ M (Q> A )l 2 ^(g, n QX + cj ) 

Q.\ 



(9) 



0, 



T 



T 



(10) 



in which n{x) = (e x — 1) _1 is the Bose distribution func- 
tion, u)q = q ■ v, and Il2(q,oj) stands for the imaginary 
part of the Fourier spectrum of the relative-electron den- 
sity correlation function defined by 

n( q ,t- = -i6{t - t')([ Pq (t), p- g (t')]) , (11) 



where p q {t) 



p q c lWor * and (...)o denotes the sta- 



electron density operator p q — ,■ e 19 ' Vj , the electron- tistical averaging over the initial density matrix p 
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The force-balance equation <j8j) describes the steady- 
state two-dimensional magnetotransport in the surface 
state of a TI. Note that the frictional forces fi and / p 
are in the opposite direction of the drift velocity v and 
their magnitudes are functions of v = \v\ only. With 
the drift velocity v = (v, 0) in the x direction, the force- 
balance equation Eq. ([5} yields a transverse resistivity 
Rxy = ~Ey/(Nev) = —B/(Ne), and a longitudinal re- 
sistivity R xx = -E x /(Nev) = -{fi + f p )/(N 2 e 2 v). The 
linear one is in the form 



1 

^£|M(Q,A)| 2 ^ 77 2 ( 9 ,^ QA ) 

2 I ^QA 



2TN 2 e 



x csch 



2T 



(12) 



III. DENSITY CORRELATION FUNCTION IN 
THE LANDAU REPRESENTATION 

For calculating the electron density correla- 
tion function II2 (<Z,w) we proceed in the Landau 



representation 
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The Landau levels of the single- 



particle Hamiltonian h = vp(n x a v — tt v (t x ) + ^g z ii-g.Bcr z 
of the relative-electron system in the absence of electric 
field are composed of a positive "+" and a negative "— " 
branch 2 ^ - — 



with e s 
0) level 



= ±y / 2ne 2 + 5 2 z 
vpsfeB and S z 



±e„ (n=l,2,....) (13) 
= — ^ffz/XB-B, and a zero (n = 



1 



-g z HBB. 



The corresponding Landau wave functions are 
1 , t 



and 



for n 




0n,fcx (y) 



n.k 



(r) = e ik ^ 





0o,kAy) 



(14) 



(15) 



(16) 



(17) 



for n — 0. Here k x is the wavevector of the system 
along x direction; 7Z n = 1 + V 2 with V n = s/2ne s / (6 Z + 
^2ne 2 + S 2 ); and 0„, fcx (y) = D n exp(-7 2 /2)ff„( 7 ) is 
the harmonic oscillator eigenfunction with H n (x) being 

the Hermite polynomial, 7 = (y — y c )/l>B = V eB(y — 
k x l%), and D n = l/(2"n!) 1 /2( eB / 7r )i/4. 



Each Landau level contains n-Q = eB/2n = l/(27rZg) 
electron states for system of unit surface area. The pos- 
itive branch e+ = e n and the n = level £0 of the above 
energy spectra are indeed quite close to those of the sur- 
face states in the bulk gap of Bi2Se3-family materials 
derived from microscopic band calculation^ 

The Landau levels are broadened due to impurity, 
phonon and electron-electron scatterings. We model the 
imaginary part of the retarded Green's function, or the 
density-of-states, of the broadened Landau level n (writ- 
ten for "+" -branch and n = levels), using a Gaussian- 
type formi2£ 



ImG„(e) 



r 



■ cxp 



2(e 



r 2 



(18) 



with a half- width J 1 of the form<^ r = [2w c /(7rr s )] 1 ^ 2 . 
Here t s is the single-particle lifetime and oj c = eBvp/sp 
is the cyclotron frequency of linear-energy-dispersion sys- 
tem with £p = 2uFV7riV being the zero-temperature 
Fermi level. Using a semi-empirical parameter a to re- 
late t s with the transport scattering time r tr = 4ar s , 
and expressing r tr with the zero-field mobility fj, at finite 
temperature^ we can write the Landau-level broadening 
as 



r = (ev F /TT)[2Ba/(Nfi)] 1/2 . 



(19) 



In the present study wc consider the case of n-doping, 
i.e. the Fermi level is high enough above the energy zero 
of the Dirac cone in the range of "+" -branch levels and 
the states of "—"-branch levels are completely filled, that 
they are irrelevant to electron transport. 

Special attention has to be paid to the n = level, 
since, depending on the direction of exchange potential 
the effective g- factor of a TI surface state, g z , can be 
positive, zero or negative . 24 ' 25 The sign and magnitude 
of the effective g-factor determines how many states of 
the zero level should be included in or excluded from the 
available states for electron occupation in the case of n- 
doping at a magnetic field, (i) If g z =0, the n = level 
center is exactly at £0 = and the system is electron-hole 
symmetric. The total number of negative energy states 
(including the states of the lower half of the n = level 
and states of the "—"-branch levels) and that of posi- 
tive energy states (including the states of the upper half 
of the n = level and states of the "+" -branch levels) 
do not change when changing magnetic field. Therefore, 
the lower-half negative energy states of this level are al- 
ways filled and the upper-half positive-energy states of 
it are available for the occupation of particles which are 
counted as electrons participating in transport in the case 
of 71-doping. (ii) For a finite positive g z > 0, the n = 
level £0 moves downward to negative energy and its dis- 
tance to the nearest "—"-branch level is 2\8 Z \ = g z [i^,B 
closer than to the nearest "+" -branch level at finite mag- 
netic field strength B. This is equivalent to the opening 
of an increasingly enlarged (with increasing B) energy 
gap between the "+" -branch states and the states of the 
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zero-level and the "—"-branch levels. The opening of a 
sufficient energy gap implies that with increasing mag- 
netic field the states in the "+" -branch levels would no 
longer shrink into the zero-level, and thus the n = 
level should be completely excluded from the conduc- 
tion band, i.e. only particles occupying the "-(-"-branch 
states are counted as electrons participating in transport 
in the case of n-doping, when the magnetic field B gets 
larger than a certain value (depending on the magnitude 
of g z ). (iii) For a finite negative g z < 0, the n = level 
£o moves upward to positive energy and an increasingly 
enlarged energy gap will be opened between the states 
of the zero- level and the "+" -branch and the states of 
"—"-branch levels, and particles occupying the n = 
level and "-(-"-branch states are electrons participating 
in transport when the magnetic field B gets larger than 
a certain value. 

As a result, the experimentally accessible sheet density 
N of electrons participating in transport is related to the 
Fermi energy e-p by the following equation valid at finite 
g z for the magnetic field B larger than a certain value: 



TV 



1 



2{irh 



J &/(<■) £lmG B (e), 



(20) 



in which /(e) = {exp[(e — £p)/T] + f } _1 is the Fermi dis- 
tribution function at temperature T and the summation 
index n goes over (1, 2, ....) for g z > 0, or (0, 1, 2, ....) for 
g z < 0. In the case of g z = 0, 



N =- 



2(tW b ) 2 



r r OO 

J def(e) ^ImG„(e)+ImGg(e) 



(21) 



valid for arbitrary magnetic field, in which ImG^e) = 
ImG o (e)0(e). 

The imaginary part of relative-electron density cor- 
relation function in the presence of a magnetic field, 
i72(g, w )i can be expressed in the Landau representation 
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2ttZ 2 ^ 

n.n 



C n , n ,(llq 2 /2)II 2 (n,n',uj), (22) 



in which the transform factor 



Cn,n' (£) 



-™ 1 Hi! 



Ti 



n-2 



Tn 2 -n 1 



(a 



(23) 



with n\ = min(n, n'), n 2 = max(n, n'), s n = 1 — 5 n _Q, 
and L™(x) being associated Laguerre polynomials. The 
Landau-representation correlation function B 2 (n, n' , w) 
in Eq. (|22[) can be constructed with the imaginary part 
of the retarded Green's function ImG„(e), or the density- 
of-statcs, of the nth Landau level a o 19 ' 20 



77 2 (n, n' , to) 



de[f(e)-f(e + oj)] 
x ImG„(e + w)ImG^(e). (24) 




FIG. 1: (Color online) The calculated resistivity R xx as a 
function of the magnetic field B having different effective g- 
factors: g z = 0, 10 and —10 for a TI surface system with 
electron sheet density N = 1.3 X 10 12 cm -2 in the cases of zero- 
magnetic-field mobility /i = 0.2m 2 /Vs (a) and fi — 0.7m 2 /Vs 
(b). Several integer- number positions of filling factor v = 
2irN/(eB) axe marked in (b). 



The summation indices n and n' in Eq. (|22p are taken 
over (1, 2, ....) for g z > 0, or (0, 1, 2, ...) for g z < 0. In the 
case of g z = 0, Eq. ([22]) still works and the summation 
indices n and n' go over (0,1,2,...) but with ImGo(e) 
replaced by ImGg(e) in Eq. 



IV. NUMERICAL RESULTS AND 
DISCUSSIONS 

Numerical calculations are performed for the magne- 
toresistivity R xx of surface state in a uniform TI Bi2Se3. 
At zero temperature the elastic scattering contributing 
to the resistivity is modeled by a Coulomb potential due 
to charged impuritiesi 28 ' 29 U(q) = n^e 2 / (2eoKq) with rii 
being the impurity density, which is determined by the 
zero-magnetic-field mobility fi. At temperatures higher 
than 50 phonon scatterings play increasingly impor- 
tant role and the dominant inelastic contribution comes 
from optical phonons. For this polar material, the scat- 
tering by optical phonons via the deformation potential 
can be neglected. Hence, we take account of inelastic 
scattering from optical phonons via Frohlich coupling: 



\M{Q)Y 



e 2 /2/(2e Q 2 )(K oo 1 — k 1 ). In the numeri- 



cal calculation we use the following parameters; 15 ' 29 ~ — 
Fermi velocity vp = 5.0 x 10 5 m/s, static dielectric con- 
stant k = 100, optical dielectric constant k x — 20, and 
phonon energy fl — 1 A meV. The broadening parameter 
is taken to be a = 3. 

Fig.[T]shows the calculated magnetoresistivity R xx ver- 
sus the magnetic field strength B for a TI surface system 
with electron sheet density N = 1.3 x 10 12 cm -2 but hav- 
ing different effective g-factors: g z = 0,10 and —10 for 
two values of zero- magnetic- field mobility /1 = 0.2m 2 /Vs 
and fi = 0.7 m 2 /Vs. representing different degree of 
Landau-level broadening. In the case without Zeeman 
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(a) 

\i = 0.2 m z /Vs/ 
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FIG. 2: (Color online) The longitudinal resistivity R xx is 
shown as a function of the magnetic field B for different values 
of zero-magnetic-field mobility: (a) /i = 0.2, (b) 0.35, (c) 0.5, 
(d) 0.65, (e) 0.8, and (f) 5m 2 /Vs. The inset of (a) illustrates 
the same for a larger magnetic-field range < B < 30 T. The 
filling factor v is plotted versus the magnetic field in (f); and 
several integer-number positions of v are also marked in (d) 
and (e). Here the surface electron density N = 1.3x 1 12 cm -2 
and the lattice temperature T — OK. 



splitting (g z = 0) the resistivity R xx exhibits almost no 
change with changing magnetic field up to 10 T, except 
the Shubnikov-dc Haas (SdH) oscillation showing up in 
the case of /i = 0.7 m 2 /Vs. This kind of magnetore- 
sistance behavior was indeed seen experimentally in the 
electron-hole symmetrical massless system of single-layer 
graphene^ In the case of a positive g-factor, g z = 10, 
the magnetoresistivity increases linearly with increasing 
magnetic field; while for a negative g-factor, g z = —10, 
the magnetoresistivity decreases linearly with increasing 
magnetic field. 

In the following we will give more detailed examina- 
tion on the linearly increasing magnctorcsistance in the 
positive g z case. 

Fig. [2] shows the calculated resistivity R xx versus the 
magnetic field strength B at lattice temperature T = K 
for system of carrier sheet density N = 1.3 x 10 12 cm -2 
and g z — 10, having different zero-field mobility (i = 
0.2, 0.35, 0.5, 0.65, 0.8 and 5 m 2 /Vs. All resistivity curves 
for mobility /i < 0.8m 2 /Vs exhibit clear linearity in the 
magnetic-field range and appear no tendency of satura- 
tion at the highest field shown in the figure. Especially, 
for the case /j, = 0.2m 2 /Vs, the linear behavior extends 
even up to the magnetic field of 30 T, as illustrated in the 
inset of Fig.[5{a). This feature contradicts the classical 
MR which saturates at sufficiently large magnetic field 

Note that here we only present the calculated R xx for 
magnetic field B larger than B c = 1 T, for which a suffi- 
cient energy gap 2\5 Z \ = g z ^B is assumed to open that 
with further increase of the magnetic field the states in 



the "+" -branch levels no longer shrink into the zero level 
and thus it should be excluded from the conduction band. 
This is of course not true for very weak magnetic field. 
When B — > the energy gap 2|<5 2 | — > 0, the situation 
becomes similar to the case of g z = 0: the whole upper 
half of the zero-level states are available to electron oc- 
cupation and we should have a flat resistivity R xx when 
changing magnetic field. With increasing B the portion 
of the zero-level states available to conduction electrons 
decreases until the magnetic field reaches B c . As a re- 
sult the resistivity R xx should exhibit a crossover from 
a flat changing at small B to positively linear increasing 
at B > B c . This is just the behavior observed in the TI 
Bi 2 Se 3 ^ 

Note that in the case of /i = 0.2m 2 /Vs, the broadened 
Landau-level widths are always larger than the neighbor- 
ing level interval: 2T > Ae„ = e n +i — e„, which requires 
fi < (4ea/N)[(y/n + 1 + \/n)/ir} 2 , even for the lowest 
Landau level n = 1, i.e. the whole Landau-level spec- 
trum is smeared. With increasing the zero-field mobility 
the magnitude of resistivity R xx decreases, and when the 
broadened Landau-level width becomes smaller than the 
neighboring level interval, 2T < Ae„, a weak SdH oscil- 
lation begin to occur around the linearly-dependent av- 
erage value of R xx at higher portion of the magnetic field 
range, as seen in Fig.[2](c), (d) and (e) for /i = 0.5,0.65 
and 0.8 m 2 /Vs. On the other hand, in the case of 
large mobility, e.g. /i = 5m 2 /Vs, where the broadened 
Landau-level widths 2T are much smaller than the neigh- 
boring level interval even for level index n as large as 30, 
the magnetoresistivity shows pronounced SdH oscillation 
and the linear-dependent behavior disappears, before the 
appearance of quantum Hall effect , 22 i 33 i 34 as shown in 

Fig.rjT). 

Abrikosov's model for the LMR requires the applied 
magnetic field large enough to reach the quantum limit 
at which all the carriers are within the lowest Landau 
levels while it is obvious that more than one Landau lev- 
els are occupied in the experimental samples in the field 
range in which the linear and non-saturating magnetore- 
sistivity was observed^ For the given electron surface 
density N = 1.3 x 10 12 cm -2 , the number of occupied 
Landau levels, or the filling factor v = 2irN/ (eB), at dif- 
ferent magnetic fields is shown in Fig.[2Jf), as well as in 
the Fig.^d) and (e), where the integer- number positions 
of v, i.e. filling up to entire v Landau levels, coincide 
with the minima of the density-of-states or the dips of 
SdH oscillation. This is in contrast with g z = case, 
where the integer number of v, which implies a filling up 
to the center position of the vth. Landau levels, locates 
at a peak of SdH oscillation, as shown in Fig.QJi. The 
observed SdH oscillations in the Bi2Se3 nanoribbon ex- 
hibiting nonsaturating surface LMR in the experiment^ 
favor the former case: a finite positive effective g z > 0. 

Next, we examine the density-dependence of the linear 
magnetoresistivity. To compare with Abrikosov's quan- 
tum magnetoresistance which suggests a R xx oc iV -2 
bchavior)^ 3 ^ we show the calculated R XX N 2 for above 



6 




N(10 12 cm" 2 ) N(10 12 cm" 2 ) 



FIG. 3: (Color online) RxxN 2 is plotted as a function of the 
surface electron density TV at magnetic field B — 3 T: (a) at 
different values of zero-field mobility fi, and (b) at different 
values of zero- field conductivity a. 



overlapped Landau level regime. 

From the above discussion, it is obvious that LMR 
shows up in the system having overlapped Landau lev- 
els and the separation of Landau levels makes the MR 
departure from the linear increase. At high tempera- 
ture, the thermal energy would smear the level separa- 
tion and phonon scatterings further broaden Landau lev- 
els. Hence, it is believed that this LMR will be robust 
against raising temperature. This is indeed the case as 
seen in Fig. [4] where we plot the calculated magnetore- 
sistivity R xx for the above system with zero-temperature 
linear mobility /i(0) = 0.6m 2 /Vs versus the magnetic 
field at different lattice temperatures. We can see that 
raising temperature to room temperature has little ef- 
fect on the linearity of MR. Due to the decreased mo- 
bility at higher temperature from phonon scattering, the 
weak SdH oscillation on the linear background tends to 
vanish. These features are in good agreement with the 
experimental report^ 
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FIG. 4: (Color online) The longitudinal resistivity of the sur- 
face state of a TI versus magnetic field B at various lattice 
temperatures. Here the zero-magnetic-field mobility at zero 
temperature is ^i(0) = 0.6 m 2 /Vs. 



LMR versus the carrier sheet density TV in Fig. [3] at fixed 
magnetic field B = 3T. The mobility is taken respec- 
tively to be (j, = 0.2,0.3,0.4,0.5 and 0.6m 2 /Vs to make 
the resistivity in the LMR regime. A clearly linear de- 
pendence of RxxN 2 on the surface density TV is seen 
in all cases, indicating that this non-saturating linear 
resistivity is almost inversely proportional to the car- 
rier density. In the figure we also show R XX N 2 versus 
TV under the condition of different given conductivity 
(7 = TVe/i = 10, 13, 16 and 20e 2 //i. In this case the half- 
width r is independent of surface density. The linear 
dependence still holds, indicating that this linear behav- 
ior is not sensitive to the modest TV-dependence of Lan- 
dau level broadening r as long as the system is in the 



V. SUMMARY 

In summary, we have studied the two-dimensional 
magnctotransport in the flat surface of a three- 
dimensional TI, which arises from the surface states with 
a wavevector-linear energy dispersion and a finite, posi- 
tive Zeeman splitting within the bulk energy gap. When 
the level broadening is comparable to or larger than 
the Landau-level separation and the conduction elec- 
trons spread over many Landau levels, a positive, dom- 
inantly linear and non-saturating magnctoresistance ap- 
pears within a quite wide range of magnetic field and 
persists up to room temperature. This remarkable LMR 
provides a possible mechanism for the recently observed 
linear magnetoresistance in topological insulator B12SC3 
nanoribbonsjiS 

In contrast to quantum Hall effect which appears in 
the case of well formed Landau levels and to Abrikosov's 
quantum magnctotransport ^£ which is limited to the ex- 
treme quantum limit that all electrons coalesce into the 
lowest Landau level, the discussed LMR is a phenomena 
of pure classical two-dimensional magnetotransport in a 
system having linear-energy-dispersion, appearing in the 
regime of overlapped Landau levels, irrespective of its 
showing up in relatively high magnetic field range. Fur- 
thermore, the present scheme deals with spatially uni- 
form case without invoking the mobility fluctuation in 
a strongly inhomogeneous system, which is required in 
the classical Parish and Littlewood model to produce a 
LMR£ 

The appearance of this significant positive-increasing 
linear magnetoresistance depends on the existence of a 
positive and sizable effective g-factor. If the Zeeman en- 
ergy splitting is quite small the resistivity R xx would 
exhibit little change with changing magnetic field. In 
the case of a negative and sizable effective g-factor the 
magnctoresistivity would decrease linearly with incrcas- 
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ing magnetic field. Therefore, the behavior of the lon- 
gitudinal resistivity versus magnetic field may provide a 
useful way for judging the direction and the size of the 
effective Zeeman energy splitting in TI surface states. 



National Basic Research Program of China (Grant 
No. 2012CB927403) and by the Program for Sci- 
cncc&Tcchnology Innovation Talents in Universities of 
Henan Province (Grant No. 2012HASTIT029). 



ACKNOWLEDGMENTS 



This work was supported by the National Science 
Foundation of China (Grant No. 11104002), the 



* cmwangsjtu@gmail.com 

1 R. Xu, A. Husmann, T. F. Rosenbaum, M. L. Saboungi, 
J. E. Enderby, and P. B. Littlewood, Nature 390, 57 
(1997). 

2 J. Hu and T. Rosenbaum, Nature Mater. 7, 697 (2008). 

3 M. P. Delmo, S. Yamamoto, S. Kasai, T. Ono, and K. 
Kobayashi, Nature 457, 1112 (2009). 

4 H. G. Johnson, S. P. Bennett, R. Barua, L. H. Lewis, and 
D. Heiman, Phys. Rev. B 82, 085202 (2010). 

5 A. L. Friedman, J. L. Tedesco, P. M. Campbell, J. C. Cul- 
bertson, E. Aifer, F. K. Perkins, R. L. Myers-Ward, J. K. 
Hite, C. R. Eddy Jr, G G. Jernigan, and D. K. Gaskill, 
Nano Lett. 10, 3962 (2010). 

6 P. L. Kapitza, Proc. R. Soc. A 123, 292 (1929). 

7 A. A. Abrikosov, Phys. Rev. B 58, 2788 (1998). 

8 A. A. Abrikosov, Europhys. Lett. 49, 789 (2000). 

9 M. M. Parish and P. B. Littlewood, Nature 426, 162 
(2003). 

10 C. L. Kane and E. J. Mele, Phys. Rev. Lett. 95, 226801 
(2005). 

11 M. Z. Hasan and C. L. Kane, Rev. Mod. Phys. 82, 3045 
(2010). 

12 X. L. Qi and S. C. Zhang, Rev. Mod. Phys. 83, 1057 (2011). 

13 Y. Xia, D. Qian, D. Hsieh, L. Wray, A. Pal, H. Lin, A. 
Bansil, D. Grauer, Y. S. Hor, R. J. Cava, and M. Z. Hasan, 
Nat. Phys. 5, 398 (2009). 

14 H. J. Zhang, C. X. Liu, X. L. Qi, X. Dai, Z. Fang and S. C. 
Zhang, Nat. Phys. 5, 438 (2009). 

15 C. X. Liu, X. L. Qi, H. J. Zhang, X. Dai, Z. Fang, and 
S. C. Zhang, Phys. Rev. B 82, 045122 (2010). 

16 H. Tang, D. Liang, R. L. J. Qiu, and X. P. A. Gao, ACS 
Nano 5, 7510 (2011). 

17 X. L. Lei and C. S. Ting, Phys. Rev. B 30, 4809 (1984); 



Phys. Rev. B 32, 1112 (1985). 

18 X. L. Lei, J. L. Birman, and C. S. Ting, J. Appl. Phys. 58, 
2270 (1985). 

19 W. Cai, X. L. Lei, and C. S. Ting, Phys. Rev. B 31, 4070 
(1985); X. L. Lei, W. Cai, and C. S. Ting, J. Phys. C: Solid 
State, 18, 4315 (1985). 

20 C. S. Ting, S. C. Ying, and J. J. Quinn, Phys. Rev. B 16, 
5394 (1977). 

21 J. W. McClure, Phys. Rev. 104, 666 (1956). 

22 Y. Zheng and T. Ando, Phys. Rev. B 65, 245420 (2002). 

23 M. Zarea and S. E. Ulloa, Phys. Rev. B 72, 085342 (2005). 

24 Z. G. Wang, Z. G. Fu, S. X. Wang, and P. Zhang, Phys. 
Rev. B 82, 085429 (2010). 

25 A. A. Taskin and Y. Ando, Phys. Rev. B 84, 035301 (2011). 

26 T. Ando, A. B. Fowler, and F. Stern, Rev. Mod. Phys. 54, 
437 (1982). 

27 X. L. Lei and S. Y. Liu, Phys. Rev. Lett. 91, 226805 (2003). 

28 C. M. Wang and F. J. Yu, Phys. Rev. B 84, 155440 (2011). 

29 D. Culcer, E. H. Hwang, T. D. Stanescu, and S. Das Sarma, 
Phys. Rev. B 82, 155457 (2010). 

30 O. Madelung, Semiconductors: Data Handbook (Springer 
Verlag, 2004). 

31 X. Zhu, L. Santos, R. Sankar, S. Chikara, C. . Howard, 
F. C. Chou, C. Chamon, and M. El-Batanouny, Phys. Rev. 
Lett. 107, 186102 (2011). 

32 Z. B. Tan, C. L. Tan, L. Ma, G. T. Liu, L. Lu, and C. L. 
Yang, Phys. Rev. B 84, 115429 (2011). 

33 Y. Zhang, Y. Tan, H. Stormer, and P. Kim, Nature 438, 
201 (2005). 

34 V. P. Gusynin and S. G. Sharapov, Phys. Rev. Lett. 95, 
146801 (2005). 

35 A. A. Abrikosov, J. Phys. A: Math. Gen. 36, 9119 (2003). 



